Inflaton potential reconstruction without slow-roll 
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We describe a method of obtaining the inflationary potential from observations which does not use 
the slow-roll approximation. Rather, the microwave anisotropy spectrum is obtained directly from 
a parametrized potential numerically, with no approximation beyond linear perturbation theory. 
This permits unbiased estimation of the parameters describing the potential, as well as providing 
the full error covariance matrix. We illustrate the typical uncertainties obtained using the Fisher 
information matrix technique, studying the A(^* potential in detail as a concrete example. 
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I. INTRODUCTION 

The determination of the initial power spectrum of per- 
turbations in the Universe is an essential step in using the 
cosmic microwave background to constrain cosmological 
parameters. Indeed, without an understanding of these 
initial perturbations, the cosmic microwave background 
in isolation says nothing about the values of parameters 
such as the Hubble parameter h and the density param- 
eter J7o. The reason is that the effect of the cosmology is 
on the dynamics of the perturbations, and a single time- 
slice, such as the perturbations at last scattering, says 
nothing about the dynamics. Usually, this problem is 
circumvented by assuming a parametrization of the ini- 
tial conditions; the cosmological parameters then enter 
via the dynamics converting these initial conditions into 
the conditions at last scattering. For example, a common 
assumption is that there is a power-law spectrum of gaus- 
sian adiabatic scalar perturbations. This is a popular 
choice because it fits current observations, and because 
it is a good approximation to the perturbations produced 
by the simplest inflation models MM , but is clearly quite 
specific since it requires four descriptive qualifiers.n 

Given a set of microwave anisotropy measurements, 
one must determine both the cosmological parameters 
and the parameters describing the initial perturbations 
simultaneously; it is not possible to do one and not the 
other. Given the description of the initial perturbations, 
one can then try to determine the model which gave rise 
to them. If the best-fit model proves to have passive per- 
turbations, meaning that the perturbations are observed 
to be entirely in their growing mode and hence are in- 
ferred to have existed since early in the Universe's evo- 
lution (in particular, to have already existed when their 



*In fact, as shortly discussed, there's actually a fifth qualifier 
as it is normally assumed that these perturbations are entirely 
in the growing mode. 



scale was considerably larger than the Hubble radius), 
then it is reasonable to believe that they arose via the in- 
fiationary mechanism I^Q , rather than being induced by 
topological defects or other causal mechanism. It would 
be further encouraging if the perturbations proved to be 
adiabatic and gaussian, because although infiation mod- 
els exist which violate those conditions, they are proper- 
ties of the simplest inflation models. 

Attempting to derive the underlying inflation model 
from observations has become known as inflaton poten- 
tial reconstruction HJ^ . Studies have focussed on models 
where inflation is driven by a single scalar field (j), moving 
in a potential V {(/)). Such models indeed give perturba- 
tions which are passive, adiabatic and gaussian, though 
they come in two types, scalar (density perturbations) 
and tensor (gravitational waves) which need not be per- 
fect power-laws. Even then they are not the most general 
class of models leading to that set of properties, because 
models where there is more than one scalar field can also 
give rise to that outcome. However, the single-field case 
appears to be the largest class of models which can be 
dealt with as a single set, where one aims to identify the 
member of the set responsible for the observations. If 
it turns out that there is no such member, then the net 
must be cast wider to include more complicated models. 
Amongst more general inflation models, there is no prob- 
lem in generating the required predictions of the spectra 
to test them one by one against the data, but there is no 
known way of looking at the observations and construct- 
ing an inflationary model which will generate the desired 
predictions. We stress again that unless a valid model for 
the initial perturbations can be found, one cannot obtain 
the cosmological parameters as there would be no way to 
compute a microwave anisotropy spectrum to compare 
with observations. 

Even under the single-field paradigm, the situation has 
not been wholely satisfactory, the reason being that the 
analytic results available for the spectra are only ap- 
proximate, having been calculated using the slow-roll 
approximation, within which results are known only to 



second-order. Therefore, even if unbiased estimates of 
the parameters describing the perturbation spectra are 
obtained from observations (such as the amphtude and 
the spectral index n), these are not translated into un- 
biased estimators of the inflation potential. This paper 
describes how this shortcoming can be overcome. 



II. DIRECT ESTIMATION OF INFLATIONARY 
PARAMETERS 

Throughout, we work within the single field paradigm. 
The traditional technique for obtaining observational 
predictions from an inflationary model is the follow- 
ing. The potential is specified as an analytical function. 
The perturbations are then computed using the slow- 
roll approximation, to give the perturbation spectra in 
parametrized form. For example, the density perturba- 
tion spectrum Sn (fc) (following the notation of Liddlc and 
Lyth Igl ) can be expanded as a Taylor series in In k as 

k 
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where k is the comoving wavenumber and k^ is an arbi- 
trary scale where the coefficients are evaluated. The slow- 
roll approximation gives the coefficients as functions of 
the potential, but only approximately, and there is a fur- 
ther approximation when the series is truncated at some 
level. Further, because the expression relating the scalar 
field value (j> to the scale k crossing the horizon is also ap- 
proximate, there is a problem of a 'drift of scales'; as we 
move from the expansion scale fc, , we begin to misiden- 
tify the (j) value corresponding to a fc value by more and 
more. 

Often all these errors are unimportant, especially for 
observations of the current quality. That for certain mod- 
els they give an error which will be significant for fu- 
ture observations has been noted by several authors [0-|| . 
One can attempt to improve things by going to the high- 
est possible order in the slow-roll expansion, which is 
unfortunately only second-order, or by taking more and 
more terms in Eq. (|l|) |lO|jll|], which does not require 
going to higher order in slow-roll. 

One does the best one can with the scalar pertur- 
bations, and also carries out a similar process for the 
tensor spectrum, which is less demanding theoretically 
as tensors are harder to detect observationally. These 
parametrized spectra are then fed into a numerical code 
(e.g. CMBFAST |12|, possibly enhanced to allow non- 
power-law spectra) to compute the microwave anisotropy 
spectrum, the Ci, as a function of those and the cosmo- 
logical parameters. 



In going in the reverse direction, starting with observa- 
tions, the standard procedure is to use the observations to 
estimate the coefficients in Eq. (|l|) , along with the cosmo- 
logical parameters. An example using current data is the 
analysis by Tegmark jl3] . As long as the parametrization 
of the perturbations was adequate, that's the job done as 
far as the cosmological parameters are concerned. How- 
ever, to obtain the inflationary potential, the approxi- 
mate slow-roll results which give the spectra in terms of 
the inflationary potential are inverted. This procedure 
is reviewed in Ref. 0], but does not yield unbiased esti- 
mates of the inflation potential. 

The numerical technology now exists to circumvent 
this problem. The key is to immediately abandon any 
attempt to make the calculations analytically. Instead, 
the perturbation spectra are obtained by numerical solu- 
tion of the relevant mode equations, which give the per- 
turbation amplitude at a particular wavenumber. The 
best formalism is that of Mukhanov jlj], and the only 
assumption is that linear perturbation theory is valid, 
which is more or less guaranteed by the fact that the 
observed (dimensionless) perturbations are order 10^^. 
The necessary ingredients to proceed are 

1. A program which can numerically solve the mode 
equations wavenumber by wavenumber. We de- 
scribed such a code in Ref. Q . This must be able 
to compute both scalar and tensor perturbations. 

2. A version of CMBFAST which is capable of taking 
arbitrary power spectra as input to produce the C'l 
curve. The output Ci is the sum of the scalar and 
tensor parts. Polarization anisotropics should be 
computed as well as temperature ones. 

We have assembled these codes into an IDL pipeline. 
The input step is to supply a parametrization of the 
potential, rather than an analytic form. In this paper 
we use the simplest version, a Taylor expansion about 
some scalar field value </>,, with the slight subtlety of 
pulling out the overall normalization as a prefactor for 
later convenience. This leads to the Cg as a function of 
the cosmological parameters and the potential parame- 
ters, i.e. C,(K, K'/K, K'Vl^, -, K no, f^B, ^a, ...) where 
primes are derivatives with respect to 0, evaluated at 0, . 

The inversion is now direct; the observed anisotropy 
spectrum is used to directly estimate the potential pa- 
rameters, which can be done in an unbiased way to gen- 
erate the best possible reconstruction. If at this stage one 
were to find that the overall best-fit model was a poor 
fit to the data, the first thing would be to try an im- 
proved parametrization of the potential and/or inclusion 
of extra cosmological parameters, and if that still fails it 
would be time to suspect that the single field paradigm 
is not correct. 

However, optimistically assuming that the best fit is 
adequate, we have our best-fit inflationary potential. But 
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FIG. 1. The traditional route from model to observables 
and back is the two-stage process at the top. The procedure 
outhned in this paper enables a direct route without approx- 
imations beyond linear perturbation theory. 



that's not all; a further advantage of this direct method 
is that it immediately gives us the covariance of the un- 
certainties on the potential parameters. For example, it 
is known that the errors on V* and VI will be highly cor- 
related. Using the old approach, these correlations would 
have to be carried through the complicated reconstruc- 
tion equations, an unpleasant enough task that in Ref. g] 
we instead used a Monte Carlo method to illustrate the 
uncertainties of a reconstruction from simulated data. 

In this approach, the consistency equation relating 
scalar and tensor perturbations (see Ref. Q for a dis- 
cussion) is automatically incorporated, being tested by 
whether there is a potential offering a satisfactory abso- 
lute fit to the data. It could of course also be tested in 
the traditional way by power spectrum fitting, but any- 
way it is unlikely that observations will be good enough 
to say anything significant. 

The two strategies are contrasted in Fig. 1. We do 
not view our new approach as replacing the traditional 
one, but rather as a next step that one would take if the 
traditional fitting proves successful, in order to obtain 
optimal results. 



variables) whatever the evolution is. As the observations 
cover a limited range of wavenumbers A;, we need only 
know V{(j)) for a limited range of (j) values about the time 
when the relevant scales cross outside the horizon during 
inflation. Our input assumption (ultimately to be tested 
against the observations) is that there is a potential T^(</>), 
and we will simply need a parametrization of it which is 
accurate enough over the desired range. 

There is one subtlety to this. The scalar wave equation 
is second order, so, in addition to the value of 0, it looks 
as if 4> is an arbitrary initial condition which needs to be 
considered as an extra parameter. However it has long 
been known that this is not the case, because scalar field 
cosmologies have an attractor behaviour whereby all ini- 
tial conditions quickly converge [Ijjl^] (indeed, during 
inflation convergence is at least exponentially fast with 
0, but even non-inflationary expansion exhibits this be- 
haviour). However, this does mean that we have to be 
sure that the simulation has run for long enough that the 
attractor is attained before the perturbations on observ- 
able scales are generated, exactly as is believed to have 
happened in the real Universe. 

As in known inflation models all observable scales cross 
outside the horizon over a very narrow range of 0, the 
simplest approach is a simple Taylor series expansion 
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where 0^ is arbitrary and can be set to zero if desired. 
We pulled out the normalization before expanding, as 
then the normalization of the C^ depends only on Vt, 
and not the other terms. In principle one could consider 
a more sophisticated expansion to try and improve the 
convergence properties such as a Fade approximant, but 
that can be assessed once actual data is available. 



B. Parameter uncertainty 



III. UNCERTAINTY AND COVARIANCE OF 
INFLATON POTENTIAL PARAMETERS 

A. Parametrizing the potential 

Although in principle a solution of the mode equation 
runs from an early initial time until the scale is well out- 
side the horizon, when considering perturbations on a 
given scale fc, the details of how the Universe expands 
are only important for a fairly brief interval around the 
time k = aH when the scale crosses outside the hori- 
zon. The reason is that while the scale is well inside 
the horizon the relevant timescales are much less than 
the expansion timescale and expansion can be neglected, 
while when scales are above the horizon the perturba- 
tions are frozen in at fixed values (in the appropriate 



Having obtained the Ci as a function of the potential 
and cosmological parameters, we can assess the likely 
accuracy with which those parameters can be found 
by a given experiment. This is carried out using the 
well-established Fisher matrix technique M,n9|, which 
amounts to taking the derivative of the Ce with respect 
to each of the parameters. The parameter uncertain- 
ties depend on the choice of 'correct' model and on the 
number of parameters allowed to vary. For illustration, 
we vary cosmological parameters about an underlying 
model with Hubble parameter h = 0.65, density param- 
eter Hq = 0.3, cosmological constant JIa = 0.7, baryon 
density fie = 0.05 and reionization optical depth t = 0.1. 

The potential we choose is the Xip'^ potential, and we 
take the epoch where the present Hubble radius equaled 
the Hubble radius during inflation as being 60 e-foldings 



before the end of inflation. Numerical solution of the 
equations of motion gives this as (/)* = 4.37 mpi. In fact 
the slow-roll approximation will work well for this poten- 
tial, and for instance can be used to show that gravita- 
tional waves should contribute about twenty percent of 
the signal at large angular scales. 

Following Zaldarriaga et al. M|, we consider a ver- 
sion of the Planck satellite which measures both tem- 
perature and polarization anisotropies, as described in 
Ref. |11[. We do not attempt to include the effects of 
foregrounds, as extensively studied recently by Tegmark 
et al. IpO], but choose to consider only one polarized 
Planck channel (in effect assuming that the polarized 
foregrounds can be removed using all the other channels) 
which, while rather approximate, yields similar results. 
The actual data, when available, will of course merit 
more sophisticated treatment. Our numbers are there- 
fore indicative only, and more importantly they would 
vary significantly if the assumed underlying model were 
changed — the quality of information available from re- 
construction depends strongly on which model (if any, of 
course) proves to be correct. 

The results are shown in Table I. The higher derivatives 
are not detected, but it is interesting to note that even 
with this very flat potential, the variation of the poten- 
tial during inflation, V', is detected at 7-sigma. However 
V" is not detected; this may seem a little surprising given 
that both the gravitational wave amplitude and the scalar 
spectral index (which depend on different combinations 
of V /V and V" /V) are in fact detectable for this poten- 
tial ||lq,0,ll]J2n] , but it turns out that the combination 
of these giving V" /V is not distinguishable from zero. In 
our approach one never needs to make the separation of 
scalars and tensors explicitly. 



C. Parameter uncertainty covariance 

The Fisher matrix technique also generates the covari- 
ances of the error estimates, and these are crucial in in- 
terpreting the observational constraints. In particular, 
the covariance matrix is essential in illustrating the re- 
constructions graphically; if correlations are ignored then 
the reconstruction deteriorates much more quickly with 
(j) than the true picture (the correlations allow for the 
fact that scale of the expansion 0* need not be the scale 
at which the observations are the most powerful). That 
our method gives the full correlation matrix of the recon- 
structed potential directly is its first key advantage over 
earlier techniques. 

To illustrate the quality of the reconstruction, we carry 
out Monte Carlo reconstructions with errors drawn ac- 
cording to the covariance matrix, and plot them against 
the true potential in Fig. 0. These reconstructions in- 
clude up to the fourth derivative; though as seen from 
Table I the higher derivatives are not detected, they can 
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TABLE I. Uncertainties for each parameter, marginaliz- 
ing over the remaining parameters. We stress that these are 
specifically for the At/)* model. These values correspond to the 
diagonal entries of the covariance matrix. There are substan- 
tial correlations between parameters, especially those describ- 
ing the potential, so the off-diagonal entries of the covariance 
matrix are significant. Of the potential parameters, only the 
magnitude and gradient of the potential are detected with any 
significance. 



still be assigned values according to their upper limit. 
We note that this potential is much less favourable for re- 
construction than ones explored previously, as it is much 
closer to the scale- invariant limit. 

The reconstructions indicate the second key advan- 
tage of the method proposed here over previous ones 
(e.g. Refs. [||J|]) — the reconstructed potentials are un- 
biased estimates of the true potential, being as likely to 
be too high as too low. In the upper panel, we see that 
the uncertainty in the overall normalization is quite large 
(ultimately due to a degeneracy in the effect of scalars 
and tensors on large angular scales). The lower panel 
shows the combination V' /V^'"^ which is primarily sen- 
sitive to density perturbations alone (indeed perfectly so 
in the slow-roll approximation), and which is much bet- 
ter determined (in particular, better than V or V' sep- 
arately). This figure allows us to see directly the range 
of (p which is constrained by the data; to highlight this 
we have indicated the values which </> takes while the mi- 
crowave anisotropies are being generated. At larger 0, 
corresponding to larger scales, the perturbations are un- 
observable, and even some way within the current hori- 
zon scale cosmic variance contributes significantly to the 
spread. Near the centre of the data the determination is 
at its best, and on short scales the information again be- 
comes poor, partly because of the dependence on all the 
cosmological parameters and partly because Silk damp- 
ing erases the perturbations as one goes beyond £ ~ 1000. 

This figure highlights once again that the information 
available from reconstruction constrains only a tiny por- 
tion of the scalar field potential. Nevertheless, the infor- 
mation available there is of good accuracy, and can be 
highly constraining in instances where theoretical mo- 
tivation suggests a potential containing few unspecified 
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FIG. 2. Twenty Monte Carlo reconstructions of the po- 
tential, compared against the true potential which is shown 
as a dashed line. The upper panel shows the potential it- 
self, and the lower one the combination V' /V'^''^ which is a 
combination coming primarily from the density perturbations 
alone. The dotted vertical lines indicate the region of the po- 
tential directly probed by the microwave background, ranging 
from the current horizon scale to the horizon scale when the 
£ — 1500 mode was generated (evaluated in the underlying 
model). The upper panel shows that the gradient is quite 
well recovered but the overall amplitude much less so, while 
the lower highlights the obvious fact that the reconstruction 
is accurate only where there is data available to constrain it. 



parameters. 

The quartic potential is an interesting test case because 
it is not far from the slow- roll limit, and this is the first 
time such a potential has been used to test reconstruc- 
tion methods. However the true strength of the method 
would be unveiled if the true model does not satisfy slow- 
roll well, despite the potential being smooth. An example 
is the potential introduced by Wang et al. 0], which was 
used to test the traditional reconstruction technique in 
Ref. [g|. It was shown in that latter paper that tradi- 
tional reconstruction could still work well, but led to a 
bias in the estimate of the potential (albeit within obser- 
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FIG. 3. Twenty Monte Carlo reconstructions of the com- 
bination V'/V^"^'^, as in the lower panel of Fig. H, but for the 
model investigated by Wang et al. |7|. 



vational errors) . Figure illustrates test reconstructions 
of this potential, using the techniques of the present pa- 
per. The uncertainties on the cosmological parameters, 
and on V and V , are very similar to those of the quar- 
tic case. However, in addition V" is detected at around 
3-sigma and the next two derivatives have uncertainties 
comparable to their values. 

IV. DISCUSSION 

It may well be that one of the simplest models of infla- 
tion is correct, and the perturbation spectra are perfectly 
satisfactorily approximated by a power-law [or at least 
some low-order truncation of Eq. (p]. If so, then cosmo- 
logical parameter estimation can proceed as described in 
previous works. However, given the intellectual and fi- 
nancial investment in pursuing cosmological parameters, 
it is vital to be aware of possible difficulties, and to an- 
alyze ways of dealing with them. We have considered 
one such possibility — that the single field paradigm is 
correct but slow-roll is not very good — and explained 
that this is readily dealt with using existing numerical 
technology. In using this technique, as with others, it is 
imperative to make an overall goodness-of-fit test to en- 
sure that the class of models being considered is capable 
of adequately explaining the data. Even if the power- 
law approximation proves valid (and certainly this is the 
method which should be tried first), one will want to 
use these techniques to ensure that estimates of the in- 
flationary potential are unbiased ones, and to obtain the 
fullest possible information about the inflaton potential 
from observations. 
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